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Abstract 

We study the relation between quantum entanglement and electron correlation in quantum 
chemistry calculations. We prove that the Hartree-Fock (HF) wave function does not violate Bell's 
inequality, thus is not entangled while the configuration interaction (CI) wave function is entangled 
since it violates Bell's inequality. Entanglement is related to electron correlation and might be used 
as an alternative measure of the electron correlation in quantum chemistry calculations. As an 
example we show the calculations of entanglement for the H2 molecule and how it is related to 
electron correlation of the system, which is the difference between the exact and the HF energies. 
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I. INTRODUCTION 



Entanglement is a quantum mechanical property that describes the nonlocal correlation 
between quantum systems which has no classical analogi>2isL££ It is a distinguishing feature 
of quantum mechanics and is regarded as a resource for many applications such as quantum 
communication and quantum computation. The basic criterion is that an entangled state 
violates Bell's inequality. A thorough understanding of entanglement has not been reached. 
The characterization of entanglement in a system of identical particles has proven particu- 
larly challenging. A number of entanglement measures for system of identical fermions have 
been proposed^^*^ 1 ^ with focus on different aspects of the phenomenon^. 

Electron correlation in a chemical system is over all electrons and is nonlocal because of 
the identical particle property of electrons. Electron correlation in quantum chemistry is 
characterized in terms of the correlation energy, which is defined as the energy difference 
between the exact non-relativistic energy of the system and the Hartree-Fock (HF) energy 
obtained in the limit that the basis set approaches completeness^. The correlation energy 
is resulted from the Coulomb repulsion and Pauli's exclusion principle. In the HF theory, 
the ground state is assumed to be described by only one Slater determinant. In dealing 
with the electron correlation effect, the configuration interaction (CI) approach postulates 
that the total wave function of the system is a linear combination of Slater determinants. 
The correlation energy is calculated as the difference between the full CI (FCI) energy and 
the HF energy using the same one-electron basis set, this is called the basis set correlation 
energy. As the one-electron basis goes to completeness, this basis set correlation energy 
approaches the exact correlation energy. 

Vedral has pointed out that interaction in general gives rise to entanglement- 4 . Here, we 
analyze the relation between the entanglement and correlation energy in molecular system. 
We show that the HF wave function, a single Slater determinant, does not violate Bell's 
inequality, thus it is unentangled. But the CI wave function violates the Bell's inequality, it 
is entangled. According to this we suggest that electron correlation is related to quantum 
entanglement. The result of the ab initio calculations for H 2 molecule is shown to support 
our suggestion. Since there is no operator in quantum chemistry for electron correlation 
energy, and measurement for quantum entanglement has been suggested, thus entanglement 
might be used as an alternative way to measure the electron correlation. 
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II. ELECTRON CORRELATION IN MOLECULAR SYSTEMS 



Electron correlation has a strong influence on many atomic, molecular—, and solid 
properties 18 of physical systems. For many-particle systems, like molecular systems, the 
state wave functions are vectors in Hilbert space of the system. Physically, these wave 
functions describe, under the Born-Oppenheimer approximation^ and the orbital approx- 
imation, how the electrons are distributed between orbitals (one-electron wave functions). 
A distribution of electrons between the orbitals, obeying the Pauli exclusion principle, de- 
termines an electron configuration. Each configuration creates a set of configuration state 
functions (CSFs). CSFs are constructed from a set of Slater determinants, the CSFs span 
the Hilbert space of the system. 

The HF method^ 1 ^ assumes that the ground state wave function is a single Slater de- 
terminant. But for most systems this is not sufficient, electron correlation needs to be con- 
sidered in order to accurately describe a molecular system. In quantum chemistry, electron 
correlation can be divided into two parts: dynamic electron correlation and non-dynamic 
electron correlation. 

If a number of electron configurations are relatively close in energy (degenerate or near- 
degenerate), the HF approximation can not provide an adequate description of the ground 
state wave function. This becomes apparent when one explores regions of Born-Oppenheimer 
potential energy surfaces(PESs) far away from the equilibrium structure. An accurate de- 
scription should consider the expansion of all these configurations. A description of the 
ground state using these degenerate or near-degenerate electron configurations comprises 
the non-dynamic electron correlation. In quantum chemistry the multiconfigurational self- 
consistent field (MCSCF) method is used to consider the non-dynamic electron correlation. 
In this method both the molecular orbitals and the configuration expansion coefficients in a 
small configuration space are optimized variationally. 

Electrons in a system tend to avoid each other because of the Coulomb repulsion, this is 
called dynamic electron correlation. Configurations created by excitations of electrons from 
the non-dynamically-correlated wave function can be considered with which to describe the 
dynamic correlation. In quantum chemistry, the CI approach is used to describe the dynamic 
electron correlation. If the wave function is expanded over all CSFs that constructed from 
all possible electron configurations, then this wave function is called the FCI wave function. 
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In the FCI approach, all the configuration coefficients are optimized variationally. 

For a given basis set, the most accurate wave function that considers both non-dynamic 
and dynamic electron correlation is a wave function in the FCI space In this wave function, 
each electron configuration will have a certain contribution to the total wave function. This 
total wave function is: 

^total = Yl A K$K, $> K = Aj[ = £ XvP\A (1) 

K iCK \i 

which is a linear combination of all the CSFs in the full CI space. Each CSF differs in how 
the electrons are distributed in the MOs, fa. The MOs are expanded in a basis of atomic 
orbitals, x M . For this wave function both the configuration expansion coefficients Ax and 
the MO expansion coefficients are variationally optimized. The optimized vector is the 
best approximation to the exact wave function of the system for a given one-electron basis 
set. 

After the FCI wave function is obtained, the optimized result can be expressed in terms of 
the reduced one particle density matrix, which gives the occupation number on each orbital 
after considering electron correlation. Thus the reduced one particle density matrix contains 
all the information about electron correlation. 

III. ENTANGLEMENT IN SYSTEMS OF IDENTICAL PARTICLES 

Entanglement in a system of identical particles is fundamentally different from that of 
the distinguishable particles, for which entanglement is invariant under local unitary trans- 
formations. In contrast, there is no local operation that only acts on one of the identical 
particles^. In the identical particle system, the single-particle basis transformation is made 
on each particle and a different set of particles is chosen in representing the many-particle 
system. Thus the entanglement property of a system of identical particles depends on the 
single-particle basis usecr^^. 

In a system of fixed number of electrons, characterization of the entanglement must 
exclude the non-factorization of the wave function due to antisymmetrization. It can be 
proved that a wave function formed through antisymmetrization of a factorized product 
state, a Slater determinant, does not violate Bell's inequality^, thus it is unentangled; while 
a CI wave function violates Bell's inequality, it is entangled. 
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For an iV-electron system, the Slater determinant is: 



|*(l,2,...,iV)> = \ Xj (l)Xk(2)... X i(N) > 

Xj(i) Xfc(i)-xKi) 

1 



Xj(2) Xk(2)... X i(2) 
Xj (N) X k(N)...Xi(N) 



(2) 



where X j s are the spin orbitals. We define an operator: 



6 = a(l)aP j {l)®a(2)bP k {2) + o(l)aPj{\) ® a(3)6P fc (3) + ... + a(l)aP,(l) ® a(2)bP k (N) + 
a(2)aP j (2) <g> a{l)bP k {\) + a(2)aPj(2) <g> a(3)6P fc (3) + ... + &{2)aP j {2) ® a(2)bP k (N) 
+... + 

ff(N)aPj(N) <8> a(l)6P fe (l) + a^aPjiN) ® <x(2)&P fe (2) + ... + 

a(N)aPj(N) ® a(JV - l)6P fc (JV - 1) (3) 

where P is the projection operator, a and 6 are any real three-dimensional unit vectors, 
the observable cxa can be refered to as a measurement of spin along the a axis, it gives the 
result +1 or —1. We can calculate the expectation value E(a,b): 



E(a,b) = <*(1,2,...,N)\0\*(1,2,...,N)> 

= < XjWXj >< Xk\°b\Xk > (4) 

in which it can be factorized as tensor product. The occurrence of a factorized product 
of two mean values implies that no choice of the unit vectors a, b, c and d can lead to a 
violation of Bell's inequality-^: 

E(a, b) + E(b, d) + E(c, d) - E(a, c)<2 (5) 
So the HF wave function is not entangled. 
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For a superposition of spin conserved Slater determinants, a spin adapted CI wave func- 
tion, it can be proved that it violates Bell's inequality. For an iV-electron system, the 
dimension of the configuration space is D, the CI wave function is: 



D 



D 



= E^lx}(l)xU2)...xj(A r )> (6) 

i 

where % represents i-th Slater determinant. We define the operator O: 

6 = *(i)TO^(i)]®*(2$©^ 

i i i i 

i i 

ff(2)S[@F}(2)] ®^(1)&[©P*(1)] + <t(2)3[©J*(2)] ® ?(3)^P'(3)] + ... + 

i i i i 

a(2)a[<$Pi(2)]®<j(2)b[<$P£(N)} 

i i 

+ ... + 

® ^(1)6[©^(1)] + Wa[0^(7V)] ® <?(2)6[©i*(2)] + ... + 

i i i i 

<W*[© ^(JV)] ® - 1)6[© JftJV - 1)] (7) 

i i 

We can calculate the expectation value E(a, b): 

D 

E(a, b) = J2 $ < Xpa\ X ) >< xi\^\xi > (8) 

i 

In general it cannot be factorized, the nonfactorization of the mean value implies that we 

— * — * 

can choose four unit vectors a, b, c, and d can lead to the violation of Bell's inequality. 

We take the ground state HF wave function of H 2 molecule as an example to show that 
this wave function does not violate Bell's inequality. The HF wave function of the ground 
state H 2 is the Slater determinant: 

An entangled state violates Bell's inequality, for example, the famous Bell state which is 
known as the spin singlet state: 



101 > -110 > 

I* > = J «» 
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where represents spin up state and 1 represents spin down state. Suppose this entangled 
pair is hold by Alice and Bob respectively, Alice performs a measurement of spin along any 
real three-dimensional unit vector a, that is, she measures the observable a. a = u x a x + 
a y a y + a z a z , where (a x , a y , a z ) are the spin Pauli matrices. The measurement will give +1 or 
— 1. Suppose Alice measures +1, then she can predict with certainty that Bob will measure 
— 1 on his qubit if he also measures spin along the a axis. Similarly if Alice measures — 1, 
she can predict with certainty that Bob will measure +1 on his qubit. Suppose Alice can 
measure the observable a. a and a.d; while Bob can measure the observable a.b and a.c. All 
these measurement will give result of +1 or —1. Alice and Bob will measure the shared pair 
of particles at the same time and choose randomly the observable to measure. Then we can 
calculate the expectation value E{ab — ac + hd + cd) and get the Bell's inequality^: 

\E(ab-ac + bd + cd) \ = \E(a,b) - E(a, c)\ + \E{b, d) + E(c, d)\ < 2 (10) 

where E(a, b) is the mean value of the product of the outcomes of two spin measurements 
along the direction a and b, this result is also known as the CHSH inequality^. However, 
for an entangled pair of particles like shown in Eq. (9), the Bell's inequality will be violated. 
For example, if they choose the following observables: 

aa = Zi, ad = Xi, ab = = , ac = = — (11) 

One can find that E{ab — ac + bd + cd) = 2\/2 > 2, violates Bell's inequality. 

Take the HF ground state wave function of H2 molecule as an example to show that this 
wave function does not violate Bell's inequality. The HF wave function of the ground state 
H 2 is the Slater determinant: 



*i(U) > 

a g (l)a(l) a' g (l)(3(l) 
a g (2)a(2) a' g (2)/3(2) 



'a 
1 

V2 



(12) 



-^[a g (l)a(l)a' g (2)P(2) - a g (2)a(2)a g (2)a(2)) 



(13) 



Here we use the restricted HF (RHF) wave function, where one spatial orbital can contain 
two spin opposite electrons, the spatial part of the two electrons are the same, but can be 



distinguished spatially, so we use a g and a to represent them. We define the operator: 

6 = a{l).aP{\) © a(2).bP'(2) + a(l).bP'(l) ® a(2).aP(2) (14) 

where P and P' are the projection operators for spatially different orbitals. It can be 
shown that the mean value E(a, b) of the product of the outcomes of two spin measurements 
along the directions a and b in two spatially different regions has the following expression: 

E{a,b) = < *|6|* > 

= l -[< a(l)\a(l)a\a(l) >< /3(2)|a(2)6|/3(2) > + 

< a(2)|5(2)a|a(2) >< /3(l)|a(l)6|/3(l) >] (15) 

because of the indistinguishable property of the electrons, the above formula can be 
written as: 

< V\d\iff >=< a\aa\a >< (3\ab\(3 > (16) 

in which it can be factorized as tensor product. The occurrence of a factorized product of 
two mean values implies that no choice of the unit vectors a, b, c and d can lead to a violation 
of Bell's inequality^. So the HF wave function is not entangled. For a superposition of Slater 
determinants, a CI wave function, it can be proved that it violates Bell's inequality. Still 
using H 2 system as an example, in the STO-3G basis, when the distance R between the two 
hydrogen atoms goes to infinity, the state can be written as: 

^(MU) > +K(U) >) = \[a g {l)cx{l)a g {2)(3{2) - a' g (l)(3(l)a 9 (2)a(2) + 

a' u {l)a{l)a u {2)(5{2) - a u (l)f3(l)a' u (2)a(2)] (17) 

For the operator: 

6 = a(l).a(P CTs (l)©P ffu (l))®a(2)i(P CT; (2)©P CT:t (2)) + 

a(l)i(P CT; (l) © P ct; (1)) © a(2).a(P CT9 (2) © P ffu (2)) (18) 

It can be shown that the mean value E(a, b) is: 
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E(a,b) 



< *|0|* > 

1 - - 

-[< a\aa\a >< (3\ab\f3 > + < a\ab\a >< (3\aa\(3 >] 



(19) 



in which the mean value can not be factorized as a tensor product, the non-factorization 
of the mean value implies that we can choose four unit vectors a, b, c, and d can lead to the 
violation of Bell's inequality^: 



So the CI wave function is entangled. It can be easily proved that any CI wave function 
is entangled. 

IV. MEASURE THE ENTANGLEMENT IN SYSTEMS OF IDENTICAL PARTI- 
CLES 

Shi^ suggested that entanglement of identical particles can be characterized by using 
the antisymmetrized basis since they are not entangled, which is equivalent to the particle 
number representation in the Fock space. This leads to the use of occupation numbers of 
different single particle basis states as the distinguishable degrees of freedom in quantifying 
entanglement of identical particles. The non-factorization due to antisymmetrization is 
naturally excluded in this representation. The occupation number entanglement in a system 
of fixed number of electrons is nothing but the situation that the state is a superposition of 
different Slater determinants. 

The particle number basis state for a fixed number of electrons is the normalized anti- 
symmetrized basis in configuration space, that is, the Slater determinants. In terms of the 
product basis \ki, k^ >, the iV-particle state is 



the coefficients q(k\, /cat) are antisymmetric, ki is the single particle state. For a fixed 
number of electrons, the normalized antisymmetrized basis can be written in terms of the 
occupation numbers of different single-particle basis states. This is the particle number 
representation, in which, 



\E(a,b) -E(a,c)\ + \E(b,d) + E(c,d)\ < 2 



(20) 




(21) 



fcl,...fcjv 
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f{ni,...,n 00 )\n 1 ,...,n 00 > 



(22) 



where nj is the occupation number of mode j, which is the spin orbital in molecules. 
The summation is subject to the constraint J2k n k = N, where N is the total number of 
electrons. 

A Slater determinant is not entangled with respect to the given single-particle basis. 
A superposition of Slater determinants is entangled in a given single-particle basis, that 
is, the spin orbital in chemical systems. A transformation from a superposition of Slater 
determinants to a single Slater determinant in another single-particle will involve operations 
on all particles in representing the state, it is a nonlocal operation. Hence, the entanglement 
in the superposition of Slater determinants is the entanglement between different single- 
particle basis states, the spin orbitals. 

The density matrix of a chemical system can be obtained from the wave function defined 
in the Fock space, from which one obtains the single particle reduced density matrices 
and the occupation numbers for each spin orbital. From the Fock space reduced density 
matrices, the bipartite entanglement between the occupation numbers of I spin orbital and 
the occupation numbers of the other spin orbitals can be calculated. 

For electrons in the molecules, the single-particle basis is the molecular spin orbital, 
which includes both spin and spatial orbital. With the spatial orbital modes as labels with 
which the particles are effectively distinguished, the entanglement can be viewed as the 
spin entanglement among the particles in different spatial orbital modes. In this way it is 
meaningful to say that the particle in a certain orbital is spin entangled with the particle in 
another orbital and this spin entanglement can be transformed into the spatial orbital mode 
entanglement. 

As we have discussed before, the entanglement in the superposition of different Slater 
determinants can be expressed as the occupation number entanglement between spin orbitals 
in a systems of fixed number of particles. Thus, entanglement is related to the electron 
correlation energy and is given by von Neumann entropy 22 



S = ~ 1^ -7T lo §2 -7T 



(23) 



k " 

where n k is the occupation number on the fc-th spatial orbital. 
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V. ENTANGLEMENT AND CORRELATION ENERGY CALCULATIONS FOR 
THE H 2 MOLECULE 



The entanglement of superposition of Slater determinants finally results in the entangle- 
ment between the molecular spin orbitals, since different Slater determinants are constructed 
by different excitations between the spin orbitals. The spin orbital entanglement character- 
izes the entanglement in the superposition of Slater Determinants. It can be characterized 
by using the von Neumann entropy of the occupation numbers of each spin orbital. This 
occupation number is between and 1. 

The entanglement between the molecular spin orbitals is defined using the eigenvalues 
Uk of the one-particle density matrix, which is the occupation of the spin orbitals. Since 
we use the restricted Hartree-Fock wave function in the calculation, two electrons share the 
same spatial orbital and only different in the spin part, the occupation of each spin orbital 
is one-half of the occupation of the spatial orbital, the occupation number is between and 
1. So the entanglement between the spin orbitals in the system can be calculated using Eq. 
(22), which characterizes the entanglement between these N/2 spin orbitals and the other 
N/2 spin orbitals, where n& is the occupation number on the fc-th spatial orbital. In fact, 
similar formula to Eq.(21) have been used in the study of the entanglement in Hubbard 
model^ and correlation entropy of the H2 molecule^. In this paper we emphasized the 
physical meaning of the formula and show that it measures the spin entanglement of the 
system. 

Electron correlation effects in the molecular system end up as a function of the one particle 
density matrix, i.e. the occupation numbers of each orbital, thus we can see that through 
one particle density matrix, entanglement and electron correlation can be connected. We 
studied the entanglement and electron correlation energy in hydrogen molecule, using two 
basis sets: the STO-3G and 6-31G** to study the relation between these two quantities. As 
shown in the figures both results support our expectation that electron correlation is related 
quantum entanglement in chemical system. 

The entanglement for the H 2 molecule as a function of R, the distance between two 
hydrogen atoms, can be evaluated using a one-particle density matrix calculated from the 
FCI wave function. Fig. 1 shows the calculation using STO-3G basis, Fig. 2 using the 
6-31G** basis set. The correlation energy is calculated as E corr = E HF — E FC i- Both the 
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entanglement and the correlation energy goes to as R — > in the STO-3G calculation. 
As R increases, the entanglement and the correlation energy increase, when R — > oo, the 
entanglement reaches its maximum 1, and the correlation energy also reaches its maximum. 
This can be explained that as R — > oo, the overlap of the two hydrogen atomic wave 
functions goes to 0, the energy difference between the two molecular orbitals in the H 2 
system, a g and a u goes to 0, the contribution to the configuration interaction from the a u 
orbital is equivalent to that of the a g orbital, so their occupations are equal, which makes 
the entanglement of the system goes to the maximum 1. So as R — > oo, the system is in a 
superposition state, 

^(k(TKU) > +k(IK(T) >) (24) 

which is entangled. This will be clear when it is written in the expanded form: 



■^(MU) > +MU) >) = iK(l)a(lK(2)/5(2)-a 9 (l)/?(lK(2)a(2) + 

<7 u (l)a(l)<7 u (2)/3(2) - <7 u (l)/3(l)<7 u (2)a(2)] 
= \(a g (l)a g (2) + r tt (lK(2))(a(l)^(2)-/?(l)a(2)) (25) 

As R decreases, the entanglement and the correlation energy decrease, the system is in 
the state: a\a g (])a g (\) > +P\cr u (l)a u (^) > where a 2 + (3 2 = 1 and a ^ (3. When a ^ (3, 
the entanglement is less than 1. When R — > 0, the entanglement reaches its minimum 0, 
and the correlation energy also reaches 0. This can be explained that as R — > 0, the overlap 
of the two hydrogen atomic wave function increases and goes to 1 when R = 0, that is the 
He atom. The energy difference between the two molecular orbitals in the H 2 system, o g 
and a u increases and goes to oo when R goes to 0, the contribution to the configuration 
interaction from the a u orbital goes to 0, so the occupation on a u orbital decreases to and 
the occupation on a g orbital increases to 2. As R = 0, for the He atom the two electrons are 
located in the Is orbital, the system is completely described by one Slater determinant, the 
entanglement of the system goes to and there is no correlation energy. We can see from the 
figure that the correlation energy changes in the same way as the entanglement changes, this 
supports our argument at the perfectly. We also studied the system using a bigger basis set 
the 6-31G** which gives us similar results, except that there are occupations on the other 
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molecular orbitals and both the entanglement and the electron correlation energy are not 
when R — > 0. We can also see that as the basis set becomes larger, the rescaled entanglement 
and the correlation energy become closer. The rescaled entanglement is simply obtained by 
multiplying it by a constant to force the entanglement at infinity to be equal to the electron 
correlation energy. 

VI. EXCHANGE CORRELATION AND ENTANGLEMENT 

Next we will discuss the relation between spin exchange correlation and quantum entan- 
glement. As R — > oo, the electron correlation energy is equal to the exchange correlation 
energy, which is the exchange integral that quantifies spin correlation in a molecular system. 
In the STO-3G calculation, the electron correlation energy is given as^: 



where S\ and e 2 are the energies of the a g and a u orbitals respectively, Jn, J 22 an d Jn 
are the Coulomb integrals for Coulomb repulsion between electrons on o g orbital, a u orbital 
and between a g orbital and a u orbital, Kyi is the exchange integral. 

When R — > oo, the energy difference between the a g and a u orbitals goes to 0, the a g and 
o u orbitals are degenerate. For these two degenerate and orthogonal orbitals, Kyi reaches 
its maximum as R — > oo, meanwhile Jn, J22 and Jn reach their minimum. The Coulomb 
integral Jn, J 22 , J12 and K 12 equal to each other. Such that when R — > 00, the correlation 
energy is equal to the exchange integral, E corr = —Kyi- So as R — > 00, the exchange integral 
measures the spin correlation and can be used to characterize the spin entanglement, which 
is related to the occupation number entanglement, which reaches the maximum as the a g 
and a u orbitals are degenerate. We expect the same behavior upon increasing the basis set 




(26) 



where, 




(27) 



to 6-31G**. 



13 



VII. SUMMARY 



The concept of electron correlation as defined in quantum chemistry calculations is useful 
but not directly observable, i.e., there is no operator in quantum mechanics that its mea- 
surement gives the correlation energy. Moreover, there are cases where the kinetic energy 
dominates the Coulomb repulsion between electrons, the electron correlation alone fails as a 
correlation measure^ 1 ^. Here, we have studied the relation between the electron correlation 
energy and the entanglement in H2 system. We show that entanglement can be used as 
an alternative measure of electron correlation in quantum chemistry calculations. The spin 
entanglement can be measured using the occupation number entanglement between spin 
orbitals. Work is underway to investigate entanglement for larger molecular systems. 
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FIG. 1: Comparison between entanglement and electron correlation energy as a function of 
internuclear distance R for the H2 molecule using the Gaussian basis set STO-3G. 
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FIG. 2: Comparison between entanglement and electron correlation energy as a function of 
internuclear distance R for the H2 molecule using the Gaussian basis set 6-31G**. 
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